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FOREWORD 


This  study  was  undertaken  as  part  of  an  investigation  of  naviga¬ 
tion  within  the  solar  system  by  optical  means.  The  objective  of  the 
investigation  is  to  evaluate  the  suitability  of  various  physical 
phenomena  as  sources  of  navigational  information  end  to  estimate  the 
accuracy  of  navigational  information  obtained  by  various  techniques. 

The  work  was  supported  under  Contract  No.  AF33 (616) -7413  by  the 
Navigational  and  Guidance  Laboratory,  Wright  Air  Development  Division, 
Air  Research  and  Development  Command,  United  States  Air  Force. 

The  present  report  is  one  of  a  series  of  three  entitled: 

Satellite  Navigation  by  Terrestrial  Occupations  of  Stars 

I.  General  Considerations  Neglecting  Atmospheric  Refraction  and 
Extinction  by  All  M.  Naqvi  (GCA  Technical  Report  62-18-A). 

II.  Considerations  Relating  to  Refraction  and  Extinction  by 
Ali  M.  Naqvi  (GCA  Technical  Report  62-21-A). 

III.  Interference  Due  to  Brightness  of  Earth's  Atmosphere  by 

Rollin  C.  Jones  and  Ali  H.  Naqvi  (GCA  Technical  Report  62-22-A). 
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The.  refraction  and  extinction  of  starlight  In  the  earth's  atmo¬ 
sphere,  including  extinction  due  to  the  clouds,  as  applied  to  satel¬ 
lite  navigation  are  discussed.  It  is  concluded  that  the  unpredictable 
variation  in  the  amount  and  height  of  clouds  make  the  observations 
of  occultations  at  the  earth's  true  edge  unsuitable.  A  method  of 
"virtual  occultation"  is  described  in  which  a  change  in  direction, 
intensity,  etc.,  of  starlight  due  to  passage  through  the  atmosphere 
yields  information  from  which  the  impact  parameter  of  the  observed 
ray  can  be  determined  and  the  satellite  orbital  elements  calculated. 
Light  extinction  is  preferred  over  refraction.  A  wave  length  band 
covering  all  or  part  of  the  region  between  0.3  and  0.55  u  is  recom¬ 
mended  for  extinction  observations.  Calculations  of  extinction  in 
the  atmosphere  are  presented  and  further  improvements  are  described. 

An  error  of  approximately  km  in  the  satellite  position  measurement 
is  anticipated  due  to  unpredictable  variations  in  atmospheric  density. 


SECTION  1 


INTRODUCTION 


In  another  report  (Naqvi,  1962),  hereafter  referred  to  as  Report 
I,  a  discussion  of  star  occultatlons  by  the  earth  for  satellite  navi¬ 
gation,  has  been  given,  neglecting  atmospheric  effects.  The  measure¬ 
ments  of  the  time  of  occultation  of  six  different  stars,  of  known 
coordinates,  is  in  principle  sufficient  for  the  determination  of  the 
satellite's  orbital  elements.  Just  before  an  occullatior*  the  starlight 
must  pass  through  the  earth's  atmosphere  where  it  undergoes  refraction, 
extinction,  et.c.  Unfortunately  the  atmosphere  is  variable;  some  of 
the  variations  are  fairly  regular  and  predictable,  but  others  are  not 
accurately  predictable.  These  variations  can  and  will  cause  signi¬ 
ficant  errors  in  the  measurements  of  the  occultation  time.  In  this 
report  we  will  investigate  the  two  most  important  atmospheric  effects 
and  arrive  at  some  suggestions  for  suitable  occultation  observations. 
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SECTION  2 


FREQUENCY,  HEIGHT  AND  TRANSMISSIVITY  OF  CLOUDS 

It  is  well  known  that  clouds  consist  of  very  small  water  droplets 
and,  at  high  altitudes,  of  ice  crystals.  They  are  divided  into  sev¬ 
eral  different  types,  whose  amount  and  height  vary  with  time  and  with 
geographical  location  in  a  very  unpredictable  manner.  The  distribu¬ 
tion  of  the  average  amount  and  height  of  various  cloud  types  for 
different  latitudes  and  seasons  is  tabulated  by  London  (1957).  Of  all 
the  different  types  the  cirrus  clouds  are  found  at  the  greatest 
heights.  Their  average  amount  varies  from  10.0"  (in  winter,  at  80°- 
90°  latitude,  and  at  an  average  height  of  6.8  km)  to  34.1%  (in  Oct¬ 
ober  at  70°-80°  latitude  and  at  an  average  height  of  5.9  km).  The 
average  height  of  the  base  of  the  clouds  varies  from  5.3  km  (in 
October,  at  80o-90°  latitude,  and  with  an  average  amount  of  25.0%)  to 
10.8  km  (in  Summer,  at  20°-30°  latitude,  and  with  an  average  amount 
of  15.6%).  However,  it  is  not  the  variability  of  cloud  characteristics 
but  their  unpredictability  which  is  liable  to  cause  large  errors  in 
occultation  observations.  It  is  generally  believed  that  clouds  occur 
up  to  heights  of  18  km. 

London  (1957)  has  also  given  some  data  on  reflectivity  and 
absorptivity  of  various  cloud  types.  Based  on  his  data  the  maximum 
and  minimum  ttansmissivities  are  given  In  Table  2.1.  The  range  in 
transmissivities  is  also  quite  large. 
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TABLE  2.1 


TRANSMISSIVITY  OF  VARIOUS  CLOUD  TYPES 


Cloud  Type 

Transmissivity  (7.) 
minimum  maximum 

Stratus 

59 

73 

Nlnbostvatus 

75 

83 

Cumulus 

77 

77 

Cumulonimbus 

80 

80 

Altostratus 

50 

56 

Cirrus 

16 

25 

5 


As  a  result  of  the  variability  of  the  above -mentioned  cloud 
characteristics,  the  exact  time  of  occultation  becomes  uncertain. 

For  example  a  star  may  disappear  behind  a  high  altitude  cloud  (such 
as  one  at  18  km)  before  the  true  occultation  by  the  earth  can  occur. 

It  is  possible  to  devise  spectroscopic  techniques  which  can  distin¬ 
guish  an  "occultation"  by  the  clouds  from  an  occultation  by  the  earth. 
However,  unless  the  cloud  height  is  known,  such  an  observation  will 
hftve  to  be  discarded,  since  the  error  in  position  determination  can 
be  anything  up  to  18  km.  We,  therefore,  conclude  that  the  method 
of  occultation  of  stars  at  the  earth's  true  edge  is  unsuitable. 


SECTION  3 


THE  METHOD  OF  VIRTUAL  OCCULT AT IONS 

In  the  occuitation  method  the  basic  observation  is  the  Measure¬ 
ment  of  time  when  the  ray  of  starlight  grazes  the  earth's  surface 
(Report  I),  i.e.  passes  at  a  known  distance  a^  from  the  earth's 
center,  where  aQ  is  the  earth's  radius.  (Departures  from  sphericity 
ot  the  earth  should  be  considered,  but  they  are  in  the  nature  of  small 
corrections,  and  will  be  ignored  here).  The  disappearance  of  the 
star  at  the  edge  of  the  earth  provides  the  required  signal. 

It  is  in  principle  possible  to  determine  the  distance  of  closest 
approach  of  the  (undeviated)  light  ray  from  the  earth's  center  (the 
so-called  impact  parameter)  from  a  change  in  one  or  more  of  the 
characteristic  properties  of  the  ray  (direction,  intensity,  etc.) 
resulting  from  the  passage  of  the  ray  through  the  earth's  atmosphere. 
In  practice  what  we  propose  to  do  is  to  determine,  either  from  theory 
or  from  experiments,  the  precise  amount  of  this  change  for  a  given 
impact  parameter  (say  corresponding  to  a  height  of  20  km  above  the 
earth's  surface,  so  as  to  avoid  the  clouds).  The  "occuitation  obser¬ 
vation"  would  consist  of  the  measurement  of  the  time  when  the  light 
ray  has  undergone  the  predetermined  change.  The  accuracy  of  the 
method  depends  upon  the  accuracy  with  which  the  changes  in  the  char¬ 
acteristics  of  the  ray  at  various  levels  in  the  atmosphere  can  be 
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determined.  We  shall  call  this  method  the  method  of  virtual  occul¬ 
ta*,  ions.  The  conditions  for  occultatlon,  and  the  occultatlon  equation 
derived  In  Report  1  are  applicable  to  this  method  as  well.  In  fact, 
in  this  method  what  we  really  do  is  to  put  to  proper  use  the  refrac¬ 
tion  and/or  the  extinction,  which  are  otherwise  sources  of  error  in 
the  occultatlon  method. 

In  this  report  we  shall  assume  spherical  symmetry  of  the  earth's 
atmosphere.  With  this  assumption  a  ray  can  be  uniquely  characterized 
by  its  impact  parameter  r^  (Figure  3.1).  All  rays  with  the  same  im¬ 
pact  parameter,  no  matter  what  their  orientation  in  space,  undergo 
the  same  change  in  their  direction,  intensity,  etc.  Note  that  the 
height  y  above  the  earth's  surface,  given  by 

y  -  rt  -  aQ  (3.1) 

will  also  be  called  impact  parameter,  when  this  can  be  done  without 
causing  any  confusion. 
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SECTION  4 
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1 
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l 
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EXTINCTION  OF  LIGHT  IN  THE  ATMOSPHERE 

4.1  INTRODUCTION 

There  are  two  main  causes  of  extinction,  absorption  and  scatter¬ 
ing.  A  third  cause,  known  as  differential  refraction,  arises  due  to 
the  fact  that  a  beam  of  starlight  diverges  on  refraction  through  the 
atmosphere,  thus  causing  a  reduction  in  flux  per  unit  area  (Menzel, 
1961,  Section  3).  However,  extinction  due  to  differential  refraction 
is  negligible  at  satellite  distance  up  to  many  thousands  of  kilometers, 
and  will  not  be  considered  here. 

4.2  ABSORPTION 

4.2.1  Introduction.  Many  atmospheric  constituents  absorb  light 
in  the  ultraviolet,  vi.slble  and  infrared  regions.  The  chief  absorbing 
molecules  are  HgO,  N^,  0^,  0^,  N^O,  C02  and  CO.  Although  the  uw sorp¬ 
tion  coef f ;'cient8  of  some  positive  and  negative  ions  are  large,  their 
number  densities  are  so  small  as  to  make  their  contribution  to  total 
absorption  negligible.  Three  important  mechanisms  for  absorption  are: 
(1)  photoionization  and  photodissociaticn,  (2)  line  absorption  by 
atom,  and  (3)  band  absorption  by  molecules,  Photoionizr lion  and 
photodissociation  play  their  most  important  role  in  the  uv  region  be¬ 
low  0.3u.  Line  absorption  contributes  very  little  to  the  total 
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absorption.  Band  absorption,  especially  the  vibration-rotation  and  the 
pure  rotation  bands,  contribute  primarily  to  absorption  at  long  wave¬ 
length  visible  and  the  infrared  regions. 

4.2.2  Selection  of  a  Suitable  Wavelength  Region.  Prediction  of 
the  absorption  of  light  in  the  earth's  atmosphere  requires  that  the 
distribution  of  absorbing  molecules  be  either  constant  or  its  var¬ 
iations  should  be  predictable.  Although  the  distribution,  as  well  as 
the  total  quantity,  of  nearly  all  atmospheric  gases  shows  some  varia¬ 
bility,  the  gases  which  show  very  large  and  only  partially  predictable 
variations  are  water  vapor,  ozone  and  carbon  dioxide. 

The  first  water  vapor  bands  of  any  consequence  appec.*  in  the  red 
region  just  below  0.7  p.  They  become  stronger  towaids  the  longer 
wavelengths  and  dominate  the  regions  0.81  -  0.99  pi  and  1.095  -  1.280  p 
and  still  longer  wavelengths.  Numerous  bands  of  carbon  dioxide  occur 
in  the  region  1.4  -  2.5  p,  but  none  for  wavelengths  shorter  than  1.4  (i. 

Ozone  has  several  absorption  bands  in  the  near  uv  and  the  visi¬ 
ble  region  (Inn  and  Tanaka,  1959).  The  Hartley  bands  system,  consist¬ 
ing  of  several  diffuse  bands,  extends  from  0.320  -  0.234  p  followed 
by  strong  continuous  absorption  to  about  0.21  p  due  to  pho^odissocia- 
tion  of  ozone.  The  absorption  coefficient  tor  wavelengths  longer  than 
0.29  p  is,  however,  very  smell.  Huggins  bands,  extending  from  0.30 
to  0.35  p,  are  very  weak,  and  the  Chappuis  bands,  from  0.45  to  0.75  p, 
are  also  very  weak,  particularly  below  0.55  p.  Thus  the  wavelength 
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region  between  approximately  0.30  and  0.55  4  is  practically  free  of 
water  vapor,  carbon  dioxide  and  ozone  absorption.  It  is  also  free 
of  absorption  due  to  other  atmospheric  gases.  Goldberg  (1954)  states 
that  "the  bands  in  the  entire  region  XX  3200  -  £700  are  so  weak  chat 
atmospheric  band  absorption  can  be  detected  only  in  the  spectrum  of 
the  setting  sun". 

Thus  from  the  viewpoint  of  predictability  of  total  absorption, 
any  region  from  0.30  to  0.55  4  appears  to  be  the  most  appropriate 
choice.  SJnce  Rayleigh  scattering  coefficients  (Section  4.3.2)  for 
shorter  wavelengths  are  much  greater  than  for  longer  wavelengths,  we, 
somewhat  arbitrarily,  have  chosen  the  region  from  0.3  to  0.5  4,  al¬ 
though  any  other  region  in  the  above-mentioned  range,  including  the 
entire  range  from  0.30  to  0.55  4,  might  be  equally  suitable. 

4.3  SCATTERING 

4.3.1  Introduction;  Mie  Scattering.  The  two  main  sources  of 
light  scattering  in  the  atmosphere  are  dust  and  aerosols  (Mie  Scatter¬ 
ing),  and  atoms  and  molecules  (Rayleigh  scattering).  Dust  and  aero¬ 
sols  are  permanent  although  minor  constituents  of  the  atmosphere.  The 
particle  sizes  range  from  about  0.01  4  to  10  4  or  even  greater. 

Junge,  Charnon  and  Mason  (1961)  and  Junge  (1961)  have  recently  made 
direct  measurements  of  particle  densltv  by  instruments  carried  in 
balloons.  Separate  instruments  were  used  for  particle  sizes  below 
0.1  4  (Aitken  nuclei  counter)  and  above  0.1  4  (collection  techniques 
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using  inertlsl  jet  impactors).  For  particle  sizes  below  0.1  ^  they 

3 

found  concentrations  of  about  1000  per  cm  between  5  and  10  km  decreas¬ 
ing  to  almost  zero  at  20  km.  For  particle  sizes  above  0.1  4  they 

found  a  maximum  concentration  of  less  than  0.1  per  cmJ  at  20  km, 

3  3 

decreasing  to  about  0.03  per  cm  at  12  km  and  to  about  0.01  per  cm 

at  25  km.  This  indicates  that  the  usual  particle  concentration  above 

20  km  altitude  is  very  small.  Rozenberg  (1960)  has  given  optical 

evidence  for  scattering  by  aerosol  near  20  km  using  the  searchlight 

technique  and  measurements  of  the  sky  brightness  in  the  sun's  almu- 

cantar.  Both  methods  show  the  existence  of  Mie  scattering  up  to  20 

km.  Evidence  for  the  presence  of  solid  particles  at  much  greater 

heights  (85  -  100  km)  is  derived  from  observations  of  noctilucent 

clouds.  Photometric  and  polarization  measurements  yLudlam,  1957; 

•2  3 

Witt,  1960)  give  a  particle  density  of  10  per  cm  and  an  optical 
depth  of  only  2  x  i0*^. 

We  conclude  the  following:  (1)  A  light  ray  which  does  not  pass 
through  the  atmosphere  at  a  height  less  than  20  km  is  suitable  not 
only  because  it  avoids  the  clouds,  but  also  because  it  avoids  essen¬ 
tially  all  the  dust  and  aerosol  particles.  (2)  For  heights  above 
20  km  and  the  spectral  region  between  say  0.3  to  0.5  p,  the  only 
significant  cause  of  extinction  is  Rayleigh  scattering. 
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4.3.2  Rayleigh  Scattering.  Ths  energy  loss  per  unit  area  due 
to  scattering,  when  a  monochromatic  and  parallel  beam  of  radiation 
travels  a  length  ds  through  a  scattering  medium  is  given  by 

dEx  *  -FxPxds  (4.1) 

where  Fj  is  the  incident  flux  and  Px  is  the  monochromatic  volume 
scattering  coefficient.  For  Rayleigh  scattering  Px  is  given  by 


32it3  (u  -  l)2 
3NX4 


(4.2) 


where  u  is  the  refractive  index  of  the  medium  and  N  is  the  number 
of  molecules  per  cm3.  When  the  medium  consists  of  several  species, 
cs  is  the  case  for  the  earth's  atmosphere,  we  have, 


P 


X 


32n3  \ 

3X4  4  N, 


(4.3) 


where  the  summation  extends  over  the  various  species.  For  the 
earth's  atmosphere  all  but  oxygen  and  nitrogen  are  minor  constituents, 
and  these  are  the  only  two  gases  which  will  be  considered.  The 
composition  of  the  atmosphere  varies  with  height  h  and  this  variation 
should  be  taken  into  account. 
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The  optical  thickness  for  vertical  incidence  at  any  height  h, 
which  we  shall  call  the  normal  optical  thickness,  is  defined  by 

m  09 

Tx<h)  -J  0x(h)dh.  (4.4) 

Deirnendjian  (1955)  has  calculated  (3.,(h)  and  T,(h)  for  various  verti- 
cal  heights  and  wavelengths.  His  results  arc  given  in  Tables  4,1 
and  4.2.  The  optical  thickness  even  at  30  km  is  very  small  and  any 
extinction  at  a  height  above  100  km  can  be  neglected. 

The  optical  thickness  of  the  earth' 3  atmosphere  for  rays  for 
which  the  Impact  parameter  is  greater  than  the  earth' „  radius,  can 
be  calculated  very  easily.  Consider  Figure  4.1,  in  which  the  devia¬ 
tion  due  to  atmospheric  refraction  of  the  ray  is  neglected.  This  is 
particularly  justified  for  rays  which  do  not  pass  through  the  lower 
atmosphere. 

Consider  a  rectangular  coordinate  system  xy,  whose  x  «yi?  is 
along  the  direction  of  the  ray  and  whose  origin  0  is  the  point  at 
which  the  light  ray  grazes  the  earth's  surface.  The  optical  thickness 
at  point  F  along  the  ray  is  defined  by 


0 1 0 0 1 1  ~  1130  0 


Figure  4.1.  Path  of  light,  ray  in  earth’s  atmosphere.  Pracoeally 
all  the  extinction  occurs  between  points  A  and  r>. 
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Equation  (4. 5) .corresponds  to  h  *  100  km. 

Due  to  symmetry  around  the  y  axis,  the  optical  thickness  values  for 
only  the  positive  values  of  x  have  been  calculated.  The  total  optical 
thickness  of  the  atmosphere  filong  a  ray,  whose  impact  parameter  is  y,  is 
given  by 

■x  (y)  '  l.}  (0,  y) .  (4.6) 

(y)  is  plotted  against  y  in  Figure  4.4. 

For  a  satellite  whose  altitude  is  more  than  100  km,  the  total  opti¬ 
cal  thickness  of  the  atmosphere  along  a  given  ray  is  given  by  the  above 
equation;  it  does  not  depend  upon  the  actual  satellite  height,  since  ex¬ 
tinction  above  a  height  of  100  km  is  negligible  The  extinction  suffered 
by  a  ray  which  traverses  the  entire  atmosphere,  is  therefore  a  function 
of  the  impact  parameter  y  only.  For  an  optically  thin  atmosphere  it  is  given 

(y)j  (4.7) 

The  percentage  extinction  for  various  impact  parameters  has  been 
calculated.  It  is  plotted  in  Figure  4.5.  Note  that  the  atmosphere  is 
not  optically  thin  for  impact  parameters  less  than  z5  km  and  16  km  for 
wavelengths  of  0.375m  and  0.520m  respectively.  If  the  :•  ,’forcement 
of  light  intensity  due  to  multiple  scattering  is  taken  into  account, 
the  extinction,  as  given  in  Figure  4  5,  would  be  found  to  be  over¬ 
estimated  However,  for  impact  parameters  up  to  20  km  the  multiple 
scattering  eifect  is  not  large  Corresponding  to  an  Impact 


Extinction  =  1  -  exp 
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F ’  4.4.  Optical  thicknesc  r^(y)  of  the  entire  a 

various  impact  parameters  and  two  wavel 
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Figure  4.5.  Extinction  due  to  Kayleigh  scattering  for  various 
impact  parameters  and  two  wavelengths. 
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parameter  of  20  km,  the  extinction  is  407.  at  X  «  0.520  n  and  887.  at 
>  ■  0.375  u.  The  total  extinction  for  any  star  of  given  spectral 
intensity  distribution  can  be  calculated  in  a  similar  manner,  lr. 
will  vary  with  the  spectral  characteristics  of  a  star,  but  its  calcu¬ 
lation  for  individual  stars  is  beyond  the  scope  of  the  present  report. 

The  values  of  extinction  for  the  above-mentioned  impact  parameter  and 
for  the  wavelength  region  0.3  to  0.5  u  are  expected  to  lie  between 
60  and  807.. 

It  should  bo  emphasized  that  the  present  investigation  is  explor¬ 
atory  in  scope.  Several  improvements  are  possible,  and  should  be 
made  to  the  calculations  presented  here.  They  are: 

1.  Departures  from  the  assumed  spherical  symmetry  of  the 
atmosphere  should  be  taken  into  account.  Due  to  lack  of  spherical 
symmetry  the  extinction  will  depend  upon  the  orientation  o£  the  uy 
with  respect  to,  say,  the  earth's  equator. 

2.  The  small  contribution  to  extinction  due  to  the  true  absorp¬ 
tion  by  atmospheric  gases  should  be  calculated. 

3.  The  model  atmosphere  used  by  Deirmendjian  is  based  upon  the 
Rand  mode]  of  the  atmosphere  (Grimminger ,  1948)  and  rocket  data  of 
Havens,  Koll  and  La  Gow  ( 1 ^52) .  Since  more  accurate  data  is  now 
available,  new  calculation  based  upon  the  most  accurate  available 
data  would  be  desirable. 

4.  Deirmendjian  considered  only  two  constituents  of  the  atmosphere, 
nitrogen  and  oxygen.  The  contribution  of  other  gases,  although  small, 
should  be  investigated. 


22 


TABLE  4.1 


THE  VOLUME  tXl INCHON  coefficient  as  a  function  of 
HEIGHT  AND  WAVELENGTH  IN  THE  EARTH'S  MOLECULAR  ATMOSPHERE. 

(after  Dei rmendlian,  1953) 


•  -1 
~  px(h)  In  cm 
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h(km) 
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TABLE  4.2 
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THE  COMPUTED  OPTICAL  THICKNESS  AS  A  FUNCTION  OF 
WAVELENGTH  AND  HEIGHT  IN  THE  EARTH'S  MOLECULAR  ATMOSPHERE. 

(after  Dei  rmcndji.au,  1955) 


‘X 


(h) 


h(km) 

X  *  375  mti 

X  *  520  mpi 

0 

0,4542 

0,1181 

2 

0,3566 

0,0927 

4 

0,2770 

0,0720 

6 

0,2127 

0,0553 

8. . 

0,1608 

0,0418 

10 

0,1196 

0,0311 

12  . 

r\  noin 

\J  y  wy  !  * 

0,02285 

16 

0,0473 

0,0123 

20 

0,0255 

0,00663 

24 

0,0138 

0,00358 

30 

0,00546 

0,00142 
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5.  Seasonal  and  latitudinal  changes  of  the  atmospheric  density 
profile  are  well-known.  Recently  Alfuth,  Dickey  and  Alsobrook  (I960, 
1961)  have  attempted  to  give  a  standard  atmospheric  density  profile 
for  various  latitudes  and  months.  For  any  accurate  predictions  of  the 
atmospheric  ext* notion  it  will  be  necessary  to  take  these  variations 
into  account . 

6  The  classical  scattering  theory  due  to  Rayleigh  breaks  down 
for  wavelengths  equal  to  or  smaller  than  the  resonance  wavelengths  of 
the  scattering  atoms  and  molecules.  This  is  llkeiy  to  change  the 
scattering  cross-sections  for  ultraviolet  radiation  (Moses,  1962). 
Further  calculations  of  quantum  mechanical  scattering  cross-sections 
are  required. 

7  Integrated  atmospheric  extinction  in  a  given  wavelength 
region  should  be  calculated  for  a  number  of  selected  bright  stars. 

All  of  the  improvements  outlined  above  can  be  incorporated  into 
future  calculations  without  any  serious  difficulty,  though  not  with¬ 
out  a  considerable  effort.  With  the  incorporation  of  these  iiupiuve- 
ments  the  major  source  of  error  remains  the  unpredictable  temporal 
and  spatial  fluctuations  of  the  atmospheric  density,  even  with  the 
seasonal  and  the  latitudinal  effects  incorporated 

In  Figures  4  6  to  4.9.  we  reproduce  four  of  the  several  curves 
by  Alfuth.  Dickey  and  Alsobrook  (I960).  All  these  curves  are  for 
one  location  (International  Falls,  Minnesota)  for  the  months  shown 
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I 

Figure  4.b.  Percentage  deviations  of  the  actual  atmospheric  density 

(from  the  standard  density  for  four  values  of  the 

cumulative  percentage  frequency  (0.157.,  25%.  7  57.  and 
99.85%) 

(Place:  International  Falla,  Minnesota;  Month,  January 

(after  Alfuth  et  al.) 
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Figure  4. 
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Percentage  deviations  of  the  actual  at<u>snneric  density 
from  the  atandard  density  for  four  vali^s  '.£  the 
cutnulntive  percentage  frequency  (0.157.  257.,  757.  and 
99.857.) 

Pla-.e:  International  Falls,  Minnesota,  Month,  April 
(after  Alfuth  et  al.) 


27 


01  DO  I  i  -Ml  'OP 


'  Figure  4.8.  Percentage  deviations  of  the  actual  atmospheric  density 

from  the  standard  density  for  four  values  of  the 

i*  cumulative  percentage  frequency  (0.157.,  257.,  757.  and 

99.857.) 

Place:  International  Falls,  Minnesota;  Month,  July 
|  (after  Alfuth  et  al.) 
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I 
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Figure  4 . 9f.  '  Percentage  deviations  of  the  actual  atmospheric  density 
.'rom  the  standard  density  for  four  values  of  cne 
cumulative  percentage  frequency  (0.15%,  25%,  75%  and 
99.85%) 

Place:  International  Falls,  Minnesota;  Month,  October 
(after  Alfuih  et  al.) 
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The  abscissae  are  the  percentage  deviation  of  the  actual  density  from 
the  standard  atmospheric  density  (expressed  as  a  polynomial)  for 
four  values  of  the  cumulative  percentage  frequency  (0.157.,  25%,  757. 
and  99.857.).  For  example,  from  Figure  4.6,  which  is  from  the  dati 
for  International  Falls,  Minnesota  fer  the  month  of  January,  the 
following  Information  can  be  obtained  for  24  km  altitude.  I'ne 
measured  density,  which  was  such  that  0.157.,  257.,  757.  and  99. 85%. of  all 
the  measurements  had  values  below  this  density,  differed  from  the 
standard  density  by  -107.,  0.47.,  5.27.  and  21.67.  respectively.  Note 
that  the  negative  sign  implies  that  the  standard  density  is  greater 
than  the  measured  density  corresponding  to  a  given  cumulative 
percentage  frequency.  A  study  of  these  figures  reveal,  that  the 
deviations  from  the  standard  density  values  are  quite  frequently  large. 

Although  a  minimum  of  six  occultation  time  measurements  are 
sufficient  for  orbit  determinations,  many  more  observations  seem 
desirable.  The  least  square  method  could  then  be  applied  and  erroneous 
observations  could  either  be  discarded  or  their  effect  on  th“  ?ol”tion 
of  the  orbital  elements  greatly  minimized.  In  any  case,  errors  of 
about  57.  in  the  density  at  altitudes  between  20  and  30  km  will  be 
very  difficult  to  avoid.  It  is  obvious  that  for  rays  whc-sc  impact 
parameter  y  is  20  km,  almost  all  the  extinction  occurs  at  the  alti¬ 
tude  region  mentioned  above.  Since  the  optical  thickness  is  roughly 
proportional  to  the  density,  this  corresponds  to  errors  of  about  57.  in 
the  calculations  of  optical  thickness.  For  a  ray  of  impact  parameter 
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20  km,  is  approximately  equal  to  2  (Figure  4.3)  and  the  error  Ar^ 
is  approximately  0.1.  Taking  the  derivative  of  Equation  (4.7),  the 
error  in  extinction  is  found  to  be, 


A  (Extinction)  *  (1  -  Extinction)  AT^ 

*  0.02 


for  f.n  extinction  of  80%.  From  Figure  (4.5),  one  can  see  that  for 
rays  with  impact  parameter  of  about  20  km,  extinction  is  a  linear 
function  of  the  impact  parameter  y.  Therefore 
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SECTION  5 


REFRACTION 


An  investigation  of  refraction  in  a  planetary  atmosphere,  with 
application  to  occupations,  has  been  described  in  another  Geophysics 
Corporation  of  America  Technical  Report  (Menzel,  1961),  For  details 
regarding  refraction,  the  reader  is  referred  to  that  report. 

As  in  the  previous  sections,  we  again  assume  a  spherically 
symmetric  atmosphere.  The  generalized  form  of  Snell's  law  is 

•ir  sin  i  *  C  ■  r^  (5.1) 

where  r  is  the  radius  vector  of  a  point  along  the  path  of  the  ray, 

P  and  i  are  the  refractive  index  of  the  medium  and  the  angle  of  inci¬ 
dence  of  the  ray  at  this  point.  The  constant  C  is  shown  in  the  Appen¬ 
dix  to  be  equal  to  the  impact  parameter  r^  of  the  ray,  provided  p  is 
unity  at  an  infinite  distance  from  the  earth.  Note  that  at  some 
point  M  for  which  i  is 


Vh 


(5.2) 
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The  differential  equation  of  a  light  ray  for  a  spherically  syrm 
metric  atmosphere  is 


(5.3: 


~  is  zero  at  the  point  "here  i  equals  zr.  It  is  snown  in  the  Appendix 
dd  L 

that  M  is  the  point  of  closest  approach,  and  that  the  path  of  the  ray 
is  symmetrical  on  both  sides  of  M  (Figure  3.1). 


If  the  polar  angle  at  M  is  denoted  by  the  deviation  suffered 


by  the  ray  in  going  from  infinity  to  M  is 


|<b„  -  — !  ,  The  total  devia- 
|  M  2) 

tion  suffered  by  a  ray  in  passing  through  the  entire  atmosphere  is 


&  »  2K>m  -  n 


(5.4) 


For  a  ray  grazing  the  earth's  surface,  and  passing  through  the 
entire  atmosphere,  the  deviation  &  is  approximately  1 0 1 3 '  (twice  the 
so-called  horizontal  refraction;  see  Chauvenet,  1891  and  Menzel .  1961). 
The  deviation  at  first  decreases  rapidly  with  the  impact  parameter, 
and  then  levels  off  as  the  refractive  index  rapidly  approaches  unity. 
The  precision  of  the  angle  measurements  using  satellite  d<  i.ne  instru¬ 
ments  is  a  few  minutes  of  arc  only,  and  it  is  expected  that  the  method 
using  extinction  of  light,  discussed  in  the  previous  section  will 
give  greater  accuracy  in  satellite  navigation. 
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APPENDIX 


SOME  GENERAL  THEOREMS  FOR  ATMOSPHERIC  REFRACTION 

The  following  theorems  hold  for  a  spherically  symmetric  plane¬ 
tary  atmosphere,  independent  of  the  manner  in  whic.i  the  refractive 
index  varies  with  the  radius  vector,  provided  it  doc.  not  become 
discontinuous  anywhere  and  approaches  unity  (for  empty  space)  &>  _n 
infinite  distance.  Although  these  theorems  are  intuitively  obvious, 
and  perhaps  are  known  to  others,  the  present  author  has  not  seen  them 
stated  in  literature.  For  wan..  of  reference,  they  arc  stated  below 
and  elementary  proofs  are  given. 

Theorem  1.  In.e  constant  in  Snell's  law  for  a  ray  of  light 

ur  sin  i  ■*  C 

equals  the  Impact  parameter  r^  of  the  ray. 

Proof:  From  Figure  A.l  it  is  easy  to  see  that 

„Lim  r  sin  i  ■  r . 
r  -)«  i 

> 

Since  u  is  assumed  to  approach  unity  as  r  ->«>, 

Ur  sin  i  -  C  ■ 

Theorem  ?.  If  at  some  point  M  on  the  path  of  the  refracted  ray 

*  l 

d? 

is  zero,  the  radial  distance  r  of  this  point  is  a  minimum.  For 

I  each  refracted  ray,  there  is  one  and  only  one  point  with  r.'ils  property. 

I 

! 
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Figure  A-l 


Proof:  Taking  derivatives  of  the  differential  equation  of  the 


light  ray,  (Equation  5.3  ),  we  get 


-r-  .  r  ILL.  .  i  +UL. 

»2  ;2  J 

*  .  ri  J  ri  . 


Since  at  point  M, 


/dr|  _ 

(d?L "  °* 


from  Equation  (5.3  )  it  follows  that 


2  2  2 
'W.  "  ri  , 


Since  r  is  always  positive,  it  represents  the  minimum,  and  never  the 
M 

maximum,  distance.  Furthermore,  since  there  is  no  real  maximum  of 
distance  (except  at  «>),  there  can  be,  at  most,  one  minimum  value  of 
the  distance.  Again  since  r  goes  from  infinity,  through  a  set  of  fin¬ 
ite  values,  to  infinity  again,  tteremust  beat  least  one  mlni.v  .j  distance. 


Theorem  3.  The  path  of  a  light  ray  is  symmetrical  on  both  sides 
df 

of  the  point  at  which  is  zero. 

Proof:  If  t  is  the  angle  of  incidence  of  a  ray  going  from  right  to 


left  at  some  point  A,  the  angle  of  incidence  when  the  same  ray  is 


£ 


going  from  left  to  right  is  n  •  i.  Therefore  the  path  of  the  ray  on 
the  two  sides  of  A  is,  in  general,  not  symmetrical.  The  necessery 
condition  for  the  symmetry  of  the  path  around  some  point,  in  a  spher¬ 
ically  symmetrical  atmosphere,  is  that  the  angle  cf  incidence  a*-  thir 

point  from  either  side  of  the  ray  Is  — .  This  condition  is  satisfied 
dr 

at  the  point  where  is  zero. 
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